Staggered pilot transmission for channel estimation and time tracking

ABSTRACT

To allow a receiving entity to derive a longer channel estimate while limiting overhead, a transmitting entity transmits a pilot on different groups of subbands in different time intervals. N subbands in the system are arranged into M non-overlapping groups. Each group includes P=N/M subbands that are uniformly distributed across the N subbands. The transmitting entity transmits the pilot on a different subband group in each time interval, and selects all M subband groups in M time intervals based on a pilot staggering pattern. The receiving entity derives (1) an initial impulse response estimate with P channel taps based on the pilot received on one subband group and (2) two longer impulse response estimates with different lengths used for data detection and time tracking. Each longer impulse response estimate may be derived by filtering initial impulse response estimates for a sufficient number of subband groups using a time-domain filter.

CROSS-REFERENCE TO RELATED APPLICATION

This application is a divisional application of U.S. patent applicationSer. No. 10/926,884, filed Aug. 25, 2004 which claims priority to U.S.Provisional Patent Application Ser. No. 60/568,324, filed May 4, 2004,which is incorporated herein by reference in its entirety.

BACKGROUND

I. Field

The present invention relates generally to data communication, and morespecifically to pilot transmission, channel estimation, and timetracking in a multi-carrier communication system.

II. Background

Orthogonal frequency division multiplexing (OFDM) is a multi-carriermodulation technique that effectively partitions the overall systembandwidth into multiple orthogonal frequency subbands. These subbandsare also referred to as tones, subcarriers, bins, and frequencychannels. With OFDM, each subband is associated with a respectivesubcarrier that may be modulated with data.

In an OFDM system, a transmitting entity processes data to obtainmodulation symbols and further performs OFDM modulation on themodulation symbols to generate OFDM symbols. The transmitting entitythen conditions and transmits the OFDM symbols via a communicationchannel. A receiving entity typically needs to obtain relativelyaccurate symbol timing in order to recover the data sent by thetransmitting entity. The receiving entity often does not know the timeat which each OFDM symbol is sent by the transmitting entity nor thepropagation delay introduced by the communication channel. The receivingentity would then need to ascertain the timing of each OFDM symbolreceived via the communication channel in order to properly perform thecomplementary OFDM demodulation on the received OFDM symbol. Thereceiving entity also needs a good estimate of the response of thecommunication channel in order to perform data detection to obtain goodestimates of the modulation symbols sent by the transmitting entity.

The transmitting entity expends system resources to support channelestimation and time tracking, and the receiving entity also consumesresources to perform these tasks. The resources used by the transmittingand receiving entities for channel estimation and time trackingrepresent overhead. Thus, it is desirable to minimize the amount ofresources expended by both the transmitting and receiving entities forthese tasks.

There is therefore a need in the art for techniques to efficientlysupport channel estimation and time tracking in an OFDM system.

SUMMARY

Techniques for performing “staggered” pilot transmission, channelestimation, and time tracking in a multi-carrier (e.g., OFDM)communication system are described herein. To allow a receiving entityto derive a longer channel estimate while limiting the amount ofresources expended for pilot transmission, a transmitting entity maytransmit a pilot on different groups of subbands in different timeintervals (e.g., different symbol periods). N subbands in the system maybe arranged into M non-overlapping groups. Each group may include P=N/Msubbands that are distributed across the N subbands. The transmittingentity may transmit the pilot on a different subband group in each timeinterval. The transmitting entity may select all M subband groups in Mtime intervals based on a pilot staggering pattern. Alternatively, thetransmitting entity may use many or most of the M subband groups indifferent time intervals, so that a substantial number of all subbandsusable for transmission in the system are used for pilot transmission indifferent time intervals. The substantial number of subbands may be, forexample, all of the usable subbands, three quarter of the usablesubbands, at least half of the usable subbands, or some othersignificant percentage of the usable subbands. The receiving entity mayderive an initial impulse response estimate with P channel taps based onthe pilot received on one subband group. The receiving entity may derivea longer impulse response estimate (with up to N channel taps) byfiltering initial impulse response estimates obtained for a sufficientnumber of different subband groups, as described below.

The receiving entity may derive two longer impulse response estimates oftwo different lengths L₁ and L₂, which may be used for datadetection/decoding and time tracking respectively, where L₁=S₁·P andL₂=S₂·P. Each longer impulse response estimate may be derived based on adifferent time-domain filter that filters S or more initial impulseresponse estimates obtained for S or more different subband groups,where S may be S₁ or S₂. For each longer impulse response estimate, thefirst P channel taps are for a “main channel”, and the remaining channeltaps are for an “excess channel”. The coefficients for each time-domainfilter may be selected based on various criteria. For example, thecoefficients for the main channel may be selected to (1) cancel theexcess channel, (2) suppress time variation in the main channel, (3)provide an unbiased estimate of the main channel, and so on. Details ofthe filtering are described below. Various aspects and embodiments ofthe invention are also described in further detail below.

BRIEF DESCRIPTION OF THE DRAWINGS

The features and nature of the present invention will become moreapparent from the detailed description set forth below when taken inconjunction with the drawings in which like reference charactersidentify correspondingly throughout and wherein:

FIG. 1 shows a block diagram of a transmitting entity and a receivingentity;

FIG. 2 shows an interlace subband structure;

FIG. 3 shows an impulse response estimate for one interlace;

FIGS. 4A through 4C show three different pilot staggering patterns;

FIG. 5 shows a process for deriving a channel estimate used for datadetection;

FIGS. 6A and 6B illustrate ambiguity in a channel impulse responseestimate due to timing uncertainty;

FIG. 7 shows a process for performing time tracking;

FIG. 8 shows a channel estimator and a time tracking unit; and

FIG. 9 shows a filter for deriving a longer impulse response estimate.

DETAILED DESCRIPTION

The word “exemplary” is used herein to mean “serving as an example,instance, or illustration.” Any embodiment or design described herein as“exemplary” is not necessarily to be construed as preferred oradvantageous over other embodiments or designs.

FIG. 1 shows a block diagram of a transmitting entity 110 and areceiving entity 150 in an OFDM system 100. Transmitting entity 110 maybe a base station or a wireless device, and receiving entity 150 mayalso be a base station or a wireless device. A base station is generallya fixed station and may also be referred to as a base transceiver system(BTS), an access point, or some other terminology. A wireless device maybe fixed or mobile and may also be referred to as a user terminal, amobile station, or some other terminology.

At transmitting entity 110, a transmit (TX) data and pilot processor 120receives different types of data (e.g., traffic/packet data andoverhead/control data) and processes (e.g., encodes, interleaves, andsymbol maps) the data to generate data symbols. As used herein, a “datasymbol” is a modulation symbol for data, a “pilot symbol” is amodulation symbol for pilot (which is data that is known a priori byboth the transmitting and receiving entities), and a modulation symbolis a complex value for a point in a signal constellation for amodulation scheme (e.g., M-PSK, M-QAM, and so on). Processor 120provides data and pilot symbols to an OFDM modulator 130.

OFDM modulator 130 multiplexes the data and pilot symbols onto theproper subbands and further performs OFDM modulation on the multiplexedsymbols to generate OFDM symbols. For each symbol period, OFDM modulator130 performs an N-point inverse fast Fourier transform (IFFT) on Nmultiplexed symbols for N total subbands and obtains a “transformed”symbol that contains N time-domain samples. Each sample is a complexvalue to be transmitted in one sample period. OFDM modulator 130 thenrepeats a portion of each transformed symbol to form an OFDM symbol thatcontains N+C samples, where C is the number of samples being repeated.The repeated portion is often called a cyclic prefix and is used tocombat inter-symbol interference (ISI) caused by frequency selectivefading. An OFDM symbol period (or simply, a symbol period) is theduration of one OFDM symbol and is equal to N+C sample periods. OFDMmodulator 130 provides a stream of OFDM symbols to a transmitter unit(TMTR) 132. Transmitter unit 132 processes (e.g., converts to analog,amplifies, filters, and frequency upconverts) the OFDM symbol stream togenerate a modulated signal, which is then transmitted from an antenna134.

At receiving entity 150, the transmitted signal from transmitting entity110 is received by an antenna 152 and provided to a receiver unit (RCVR)154. Receiver unit 154 processes (e.g., filters, amplifies, frequencydownconverts, and digitizes) the received signal and provides a streamof input samples. An OFDM demodulator (Demod) 160 performs OFDMdemodulation on the input samples and provides received data and pilotsymbols. A detector 170 performs data detection (e.g., equalization ormatched filtering) on the received data symbols with a channel estimatefrom a channel estimator 172 and provides detected data symbols, whichare estimates of the data symbols sent by transmitting entity 110. Areceive (RX) data processor 180 processes (e.g., symbol demaps,deinterleaves, and decodes) the detected data symbols and providesdecoded data. In general, the processing by OFDM demodulator 160 and RXdata processor 180 is complementary to the processing by OFDM modulator130 and TX data and pilot processor 120, respectively, at transmittingentity 110.

Channel estimator 172 derives impulse response estimates based on thereceived pilot symbols from OFDM demodulator 160 and further derivesfrequency response estimates used by detector 170. A synchronizationunit 162 performs time tracking and determines symbol timing based onthe impulse response estimates from channel estimator 172. OFDMdemodulator 160 performs OFDM demodulation based on the symbol timingfrom unit 162.

Controllers 140 and 190 direct operation at transmitting entity 110 andreceiving entity 150, respectively. Memory units 142 and 192 providestorage for program codes and data used by controllers 140 and 190,respectively.

Data and pilot may be transmitted in various manners in system 100. Forexample, data and pilot may be transmitted (1) simultaneously in thesame symbol period using frequency division multiplexing (FDM), (2)sequentially in different symbol periods using time divisionmultiplexing (TDM), or (3) using a combination of FDM and TDM. The Ntotal subbands may also be used for data and pilot transmission invarious manners. An exemplary data/pilot transmission scheme isdescribed below.

FIG. 2 shows an interlace subband structure 200 that may be used fordata and pilot transmission in system 100. System 100 has an overallsystem bandwidth of BW MHz, which is partitioned into N orthogonalfrequency subbands using OFDM. Each subband has a bandwidth of BW/N MHz.Of the N total subbands, only U subbands may be used for data and pilottransmission, where U≦N, and the remaining G=N−U subbands may be unusedand serve as guard subbands. As a specific example, system 100 mayutilize an OFDM structure with N=4096 total subbands, U=4000 usablesubbands, and G=96 guard subbands. For simplicity, the followingdescription assumes that all N subbands may be used for data and pilottransmission. These N subbands are assigned indices of k=1 . . . N.

The N total subbands may be arranged into M “interlaces” or disjointsubband groups. The M interlaces are disjoint or non-overlapping in thateach of the N total subbands belongs in only one interlace. Eachinterlace contains P subbands, where P·M=N. The M interlaces are givenindices of m=1 . . . M, and the P subbands in each interlace are givenindices of p=1 . . . P.

The P subbands for each interlace may be uniformly distributed acrossthe N total subbands such that consecutive subbands in the interlace arespaced apart by M subbands. Each interlace m, for m=1 . . . M, mayinclude P subbands with the following k indices:(p−1)·M+m, for p=1 . . . P.  Eq (1)As shown in FIG. 2, interlace 1 contains subbands with indices k=1, M+1,2M+1, and so on, interlace 2 contains subbands with indices k=2, M+2,2M+2, and so on, and interlace M contains subbands with indices k=M, 2M,3M, and so on. The P subbands in each interlace are thus interlaced withthe P subbands in each of the other M−1 interlaces. Each interlace isfurther associated with a staggering phase m, which is equal to theindex k of the first subband in the interlace.

In general, system 100 may utilize any OFDM structure with any number oftotal, usable, and guard subbands. Any number of interlaces may also beformed. Each interlace may contain any number of subbands and any one ofthe N total subbands. The interlaces may contain the same or differentnumbers of subbands. For clarity, the following description is for theinterlace subband structure shown in FIG. 2 with M interlaces and eachinterlace containing P uniformly distributed subbands. This interlacesubband structure provides several advantages. First, frequencydiversity is achieved since each interlace contains subbands taken fromacross the entire system bandwidth. Second, a receiving entity mayrecover data/pilot symbols sent on a given interlace by performing apartial P-point FFT instead of a full N-point FFT, which can simplifythe processing at the receiving entity.

A communication channel between transmitting entity 110 and receivingentity 150 in OFDM system 100 may be characterized by either atime-domain channel impulse response or a corresponding frequency-domainchannel frequency response. As used herein, and which is consistent withconventional terminology, a “channel impulse response” or “impulseresponse” is a time-domain response of the channel, and a “channelfrequency response” or “frequency response” is a frequency-domainresponse of the channel. In a sampled-data system, the channel frequencyresponse is the discrete Fourier transform (DFT) of the channel impulseresponse. This relationship may be expressed in matrix form, as follows:H _(N×1) =W _(N×N) ·h _(N×1) and h _(N×1) =W _(N×N) ^(H) ·H _(N×1),  Eq(2)

-   -   where h _(N×1) is an N×1 vector for the impulse response of the        communication channel;    -   H _(N×1) is an N×1 vector for the frequency response of the        communication channel;    -   W _(N×N) is an N×N Fourier matrix; and    -   “H” denotes a conjugate transpose.        The Fourier matrix W _(N×N) is defined such that the (l,n)-th        entry, W_(N) ^(l,n), is given as:

$\begin{matrix}{{W_{N}^{l,n} = {\mathbb{e}}^{{- j}\; 2\;\pi\frac{{({l - 1})}{({n - 1})}}{N}}},{{{for}\mspace{14mu} l} = {{1\mspace{14mu}\ldots\mspace{14mu} N\mspace{14mu}{and}\mspace{14mu} n} = {1\mspace{14mu}\ldots\mspace{14mu} N}}}} & {{Eq}\mspace{14mu}(3)}\end{matrix}$where l is a row index and n is a column index.

The channel impulse response h _(N×1) is composed of N channel taps,with each channel tap h_(l) being defined by a zero or non-zero complexgain value at a specific tap delay l. The channel frequency response H_(N×1) is composed of N channel gains for the N total subbands, witheach channel gain H_(k) being a complex gain value for a specificsubband k.

If pilot symbols are transmitted on the P subbands in interlace m, thenthe received pilot symbols for this interlace may be expressed as:Y _(m) =H _(m) ∘X _(m) +N _(m),  Eq (4)where

-   -   X _(m) is a P×1 vector with P pilot symbols sent on the P        subbands in interlace m;    -   Y _(m) is a P×1 vector with P received pilot symbols obtained by        the receiving entity for the P subbands in interlace m;    -   H _(m) is a P×1 vector for the actual channel frequency response        for interlace m;    -   N_(m) is a P×1 noise vector for the P subbands in interlace m;        and    -   “∘” denotes the Hadamard product, which is an element-wise        product, where the i-th element of Y _(m) is the product of the        i-th elements of X _(m) and H _(m).        The vector H _(m) contains only P entries of the vector H _(N×1)        for the P subbands in interlace m. For simplicity, the noise N        _(m) is assumed to be additive white Gaussian noise (AWGN) with        zero mean and a variance of σ².

An initial frequency response estimate may be obtained for interlace m,as follows:Ĥ _(m) =Y _(m) /X _(m) =H _(m) +N _(m) /X _(m),  Eq (5)

-   -   where        -   Y _(m)/X _(m)=[y_(m,1)/p_(m,1) . . . y_(m,P)/p_(m,P)], and            y_(m,i) and p_(m,i) are respectively the received and            transmitted pilot symbols for the i-th subband in interlace            m; and        -   Ĥ _(m) is a P×1 vector for the initial frequency response            estimate for interlace m.            Ĥ _(m) contains P channel gain estimates for the P subbands            in interlace m, which may be obtained based on P            element-wise ratios of the received pilot symbols to the            transmitted pilot symbols, as shown in equation (5). If            interlace m contains unused subbands with no received pilot            symbols, then extrapolation, interpolation, and/or some            other technique may be used to estimate the channel gains            for these unused subbands.

A P-tap impulse response estimate using interlace m may be obtained byperforming a P-point IFFT on the initial frequency response estimate Ĥ_(m), as follows:ĥ _(m) =W _(m) ·W _(P×P) ^(H) ·Ĥ _(m),  Eq (6)where ĥ _(m) is a P×1 vector for the impulse response estimate forinterlace m;W _(P×P) is a P×P Fourier matrix with elements defined as shown inequation (3); andW _(P×P) is a P×P diagonal matrix containing W_(N) ^(−m,p) for the p-thdiagonal element, for p=1 . . . P, and zeros elsewhere, where

$W_{N}^{{- m},p} = {{\mathbb{e}}^{j\; 2\;\pi\frac{{({m - 1})}{({p - 1})}}{N}}.}$The channel component in the P elements of vector W _(P×P) ^(H)·Ĥ _(m)contains a phase ramp which may be expressed as: h_(m,p)=h_(p)·W_(N)^(m,p), for p=1 . . . P. The slope of the phase ramp is determined bythe staggering phase m of interlace m. The phase ramp may be removed bymultiplying each element of W _(P×P) ^(H)·Ĥ _(m) with W_(N) ^(−m,p) toobtain a corresponding element of ĥ _(m). The P elements of ĥ _(m) maybe expressed as: h_(p)=h_(m,p)·W_(N) ^(−m,p), for p=1 . . . P.

ĥ _(m) contains P channel taps and is obtained based on Ĥ _(m), whichcontains P channel gain estimates for the P subbands in interlace m.Since the actual channel impulse response h _(N×1) is composed of Nchannel taps, the initial impulse response estimate ĥ _(m) isundersampled in the frequency domain by the P subbands in interlace m.This undersampling in the frequency domain causes aliasing of thechannel impulse response h _(N×1), in the time domain. The initialimpulse response estimate ĥ _(m) may be expressed as:

$\begin{matrix}{{{\underset{\_}{\hat{h}}}_{m} = {{\sum\limits_{s = 1}^{M}{{\underset{\_}{h}}_{s} \cdot W_{M}^{s,m}}} + \underset{\_}{n}}},} & {{Eq}\mspace{14mu}(7)}\end{matrix}$

-   -   where h _(N×1)=[h ₁ ^(T) h ₂ ^(T) . . . h _(M) ^(T)]T is the        full-length actual channel impulse response;    -   h _(s), for s=1 . . . M, is a P×1 vector containing P channel        taps in h _(N×1) with tap indices of (s−1)·P+1 through s·P;    -   n is a P×1 vector of noise for the initial impulse response        estimate ĥ _(m);

${W_{M}^{s,m} = {\mathbb{e}}^{{- j}\; 2\;\pi\frac{{({s - 1})}{({m - 1})}}{M}}};{and}$“T” denotes a transpose.An “aliasing pattern” corresponding to staggering phase m may be definedas {W_(M) ^(s,m)}, for s=1 . . . M, and include the coefficients usedfor equation (7). The full-length actual channel impulse response h_(N×1) is composed of M segments. Each segment s contains P consecutivechannel taps in h _(N×1) and is represented by a vector h _(s). Equation(7) indicates that the M segments alias and combine when undersampled inthe frequency domain, and the combining coefficients are given by thealiasing pattern.

FIG. 3 shows an impulse response estimate 300 obtained based on pilotsymbols received on P subbands in one interlace. The full-length channelimpulse response h _(N×1) is composed of N channel taps with indices of1 through N. The first P channel taps in h _(N×1) are contained in h₁and may be referred to as the main channel. The remaining N−P channeltaps in h _(N×1) are contained in h ₂ through h _(M) and may be referredto as the excess channel. The excess channel taps alias whenundersampled in the frequency domain. The aliasing results in the excesschannel taps at indices of P+l, 2P+l, . . . , and (M−1)·P+l allappearing at tap index l, for l=1 . . . P. The P channel taps in ĥ _(m)thus contain P main channel taps as well as N−P excess channel taps.Each aliased excess channel tap causes error in the estimation of thecorresponding main channel tap.

A longer impulse response estimate with more than P channel taps may beobtained by transmitting pilot symbols on multiple interlaces. Oneinterlace may be used for pilot transmission in each symbol period, anddifferent interlaces may be used for pilot transmission in differentsymbol periods. The use of multiple interlaces for pilot transmissionallows the receiving entity to obtain a longer channel estimate, whichmay improve performance. By using all M interlaces for pilottransmission, it is possible to estimate the entire full-length channelimpulse response with N channel taps.

The specific interlace to use for pilot transmission in each OFDM symbolperiod may be determined by a pilot staggering pattern. Variousstaggering patterns may be used for pilot transmission. In anembodiment, a staggering pattern may select one interlace for pilottransmission in each symbol period based on the following:m _(t)=[(m _(t−1)−1+Δm)mod M]+1, with (Δm,M)=1,  Eq (8)where

-   -   t is an index for symbol period;    -   Δm is the difference between interlace indices for two        consecutive symbol periods;    -   m_(t) is the interlace to use for pilot transmission in symbol        period t; and    -   (x,y)=1 means that x and y are relatively prime (i.e., the        greatest common divisor for both x and y is one).        The −1 and +1 in equation (8) account for an interlace index        numbering scheme that starts with ‘1’ instead of ‘0’. The        interlace used for the first symbol period is m₁, where m₁ε{1 .        . . M}. Different “complete” staggering patterns may be formed        with different values of Δm. A complete staggering pattern is        one that selects all M interlaces for pilot transmission, e.g.,        in M symbol periods. As an example, with Δm=1, the M interlaces        are selected in sequential order, and the staggering pattern may        be expressed as {1, 2, 3, . . . , M}. For the case with M=8,        values of 1, 3, 5, and 7 may be used for Δm to obtain different        complete staggering patterns. Of these four values, 7 is        equivalent to 1 (in terms of performance) since Δm=1 is an        increment of one and Δm=7 is a decrement of one, and 5 is        equivalent to 3 for the same reason.

FIG. 4A shows a complete staggering pattern 400 that may be used forpilot transmission. The vertical axis represents interlace indices, andthe horizontal axis represents time. For this example, M=8 and oneinterlace is used for pilot transmission in each symbol period.Staggering pattern 400 is generated with Δm=1 in equation (8), and thecomplete staggering pattern may be expressed as {1, 2, 3, 4, 5, 6, 7,8}. The pilot is thus transmitted on interlace 1 in symbol period 1,then interlace 2 in symbol period 2, and so on, then interlace 8 insymbol period 8, then back to interlace 1 in symbol period 9, and so on.All eight interlaces are used for pilot transmission in each 8-symbolperiod duration.

FIG. 4B shows a complete staggering pattern 410 that may also be usedfor pilot transmission. Again, M=8 and one interlace is used for pilottransmission in each symbol period. Staggering pattern 410 is generatedwith Δm=3 in equation (8), and the complete staggering pattern may beexpressed as {1, 4, 7, 2, 5, 8, 3, 6}. The pilot is thus sent oninterlace 1 in symbol period 1, then interlace 4 in symbol period 2,then interlace 7 in symbol period 3, and so on. Again, the pilot istransmitted on all eight interlaces in each 8-symbol period duration.Over three symbol periods, staggering pattern 410 selects interlaceswith relative offsets of {1, 4, 7} while staggering pattern 400 selectsinterlaces with relative offsets of {1, 2, 3}. Staggering pattern 410 isthus more “spread out” than staggering pattern 400 and may providebetter performance.

FIG. 4C shows a complete staggering pattern 420 that does not satisfyequation (8) but may also be used for pilot transmission. This completestaggering pattern may be expressed as {1, 5, 2, 6, 3, 7, 4, 8}. Thepilot is transmitted on all eight interlaces in each 8-symbol periodduration.

In general, the pilot may be transmitted on any number of interlaces andon any one of the M interlaces in each symbol period. The particularinterlace to use for pilot transmission in each symbol period may beselected based on any staggering pattern, three of which are shown inFIGS. 4A through 4C. The pilot may be transmitted on all M interlacesusing a complete staggering pattern or on a subset of the M interlacesusing a “partial” staggering pattern.

A longer impulse response estimate {tilde over (h)} _(L×1)(t) with Lchannel taps, where P<L≦N, may be obtained by filtering multiple P-tapinitial impulse response estimates obtained for multiple interlaces.This time-domain filtering may be performed, e.g., with a finite impulseresponse (FIR) filter, as follows:

$\begin{matrix}{{{{\overset{\sim}{h}}_{s,l}(t)} = {\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{\alpha_{s,l}(i)} \cdot {{\hat{h}}_{l}\left( {t - i} \right)}}}},{{{for}\mspace{14mu} l} = {{1\mspace{14mu}\ldots\mspace{14mu} P\mspace{14mu}{and}\mspace{14mu} s} = {1\mspace{14mu}\ldots\mspace{14mu} S}}},} & {{Eq}\mspace{14mu}(9)}\end{matrix}$

-   -   where ĥ(t)=[ĥ₁(t)ĥ₂(t) . . . ĥ_(P)(t)]^(T) is an initial impulse        response estimate obtained for symbol period t based on a pilot        received on interlace m_(t);    -   {tilde over (h)} _(s)(t)=[{tilde over (h)}_(s,1)(t){tilde over        (h)}_(s,2)(t) . . . {tilde over (h)}_(s,P)(t)] is a P×1 vector        that is an estimate of the channel impulse response h _(s)(t)        for segment s in symbol period t;    -   α_(s,l)(i) is a coefficient for the i-th filter tap used to        derive the l-th channel tap in segment s;    -   N_(f) is the number of non-causal taps for the time-domain        filter; and    -   N_(b) is the number of causal taps for the time-domain filter.        The L-tap impulse response estimate {tilde over (h)} _(L×1)(t)        is composed of S segments and may be given as: {tilde over (h)}        _(L×1)(t)=[{tilde over (h)} ₁ ^(T)(t){tilde over (h)} ₂ ^(T)(t)        . . . ĥ _(S) ^(T)(t)]^(T), where S>1 and L=S·P. Each segment s,        for s=1 . . . S, contains P channel taps that are included in        the vector {tilde over (h)} _(s)(t). {tilde over (h)} _(s)(t) is        an estimate of h _(s)(t), which is the actual channel impulse        response for segment s.

Equation (9) indicates that the P channel taps for each segment s may beobtained by filtering N_(f)+N_(b) initial impulse response estimatesĥ(t+N_(f)) through ĥ(t−N_(b)+1), which may be obtained over N_(f)+N_(b)symbol periods for N_(f)+N_(b) different interlaces. The initial impulseresponse estimate ĥ(t) for the current symbol period t is aligned atfilter tap i=0. Equation (9) also indicates that each channel tap {tildeover (h)}_(s,l)(t) in {tilde over (h)} _(L×1)(t) may be obtained bymultiplying N_(f)+N_(b) channel taps ĥ_(l)(t−N_(b)+1) throughĥ_(l)(t+N_(f)) with N_(f)+N_(b) coefficients α_(s,l)(N_(b)−1) throughα_(s,l)(−N_(f)), respectively, and combining the N_(f)+N_(b) resultantproducts.

In general, the coefficients for each channel tap {tilde over(h)}_(s,l)(t) of each segment s may be selected separately. Furthermore,N_(f) and N_(b) may be selected for each channel tap of each segment s.For simplicity, one set of N_(f)+N_(b) coefficients may be used for allP channel taps in each segment, and S sets of coefficients may bedefined for the S segments of {tilde over (h)} _(L×1)(t) In this case,the coefficients {α_(s)(i)} for each segment s are not a function ofchannel tap index.

The time-domain filtering may also be performed using other types offilter, such as an infinite impulse response (IIR) filter. Thetime-domain filtering may also be performed using a causal filter (withN_(f)=0 and N_(b)≧1), a non-causal filter (with N_(f)≧1), or a filterwith both causal and non-causal taps. For clarity, the followingdescription is for the time-domain filter shown in equation (9).

1. Channel Impulse Response Estimate of Length 2P

To obtain a longer impulse response estimate {tilde over (h)} _(2P×1)(t)with L=2P channel taps, the initial impulse response estimate ĥ(t)obtained in symbol period t for one interlace may be expressed as:ĥ (t)= h ₁(t)+ h ₂(t)·W _(M) ^(m) ^(l) +n (t),  Eq (10)where

$W_{M}^{m_{t}} = {{\mathbb{e}}^{{- j}\; 2\;\pi\frac{({m_{t} - 1})}{M}}.}$Equation (10) is derived based on equation (7) and assumes that segments3 through M contain channel taps with zero magnitude. The vector h ₁(t)contains the first P channel taps in h _(N×1)(t) for the main channel.The vector h ₂ (t) contains the next P channel taps in h _(N×1)(t) forthe excess channel.

The coefficients for the time-domain filter for the main channelestimate {tilde over (h)} ₁(t) may be selected based on variousconstraints such as:

Cancel excess channel:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{\alpha_{1}(i)} \cdot W_{M}^{m_{t - i}}}} = 0},} & {{Eq}\mspace{14mu}\left( {11a} \right)}\end{matrix}$Suppress time variation:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{\alpha_{1}(i)} \cdot i}} = 0},} & {{Eq}\mspace{14mu}\left( {11b} \right)}\end{matrix}$Provide unbiased estimate:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{\alpha_{1}(i)}} = 1},} & {{Eq}\mspace{14mu}\left( {11c} \right)}\end{matrix}$andMinimize noise variance:

$\begin{matrix}{{\min\left( {\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{\alpha_{1}^{2}(i)}} \right)},} & {{Eq}\mspace{14mu}\left( {11d} \right)}\end{matrix}$where m_(t−i) is the interlace used for pilot transmission in symbolperiod t−i, which corresponds to the i-th filter tap. An unbiasedestimate is one for which the mean of the estimate (over noise) is equalto the perfect channel value.

Equation (11b) cancels the linear component of the channel variationover the N_(f)+N_(b) symbol periods, which would be the dominantcomponent at low speeds and/or small N_(f)+N_(b). The first constraintin equation (11a) cancels the contribution from the excess channel h ₂(t), so that {tilde over (h)} ₁(t) contains mostly components from themain channel {tilde over (h)}_(1,l)(t). The second constraint inequation (11b) suppresses time variation in the main channel h ₁(t)across the N_(f)+N_(b) symbol periods. The third constraint in equation(11c) provides an unbiased estimate of h ₁(t), so that the expectedmagnitude of {tilde over (h)}_(1,l)(t) is equal to h_(1,l)(t). Thefourth constraint in equation (11d) minimizes the noise variance in themain channel estimate {tilde over (h)} ₁(t). The number of taps(N_(f)+N_(b)) for the time-domain filter determines (1) the number ofdegrees of freedom for selecting the coefficients and (2) the number ofconstraints that may be applied in selecting the coefficients.

The coefficients for the time-domain filter for the excess channel{tilde over (h)} ₂(t) may be selected based on the various constraintssuch as:

Cancel main channel:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{\alpha_{2}(i)}} = 0},} & {{Eq}\mspace{14mu}\left( {12a} \right)}\end{matrix}$Suppress time variation of main channel:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{\alpha_{2}(i)} \cdot i}} = 0},} & {{Eq}\mspace{14mu}\left( {12b} \right)}\end{matrix}$Suppress time variation of excess channel:

$\begin{matrix}{{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{{\alpha_{2}(i)} \cdot {\mathbb{i}}}\; W_{M}^{m_{t - i}}}} = 0},} & {{Eq}\mspace{14mu}\left( {12c} \right)}\end{matrix}$andProvide unbiased estimate:

$\begin{matrix}{{\sum\limits_{i = {- N_{f}}}^{N_{b} - 1}{{\alpha_{2}(i)} \cdot W_{M}^{m_{t - i}}}} = 1.} & {{Eq}\mspace{14mu}\left( {12d} \right)}\end{matrix}$The first constraint in equation (12a) cancels the contribution from themain channel h ₁(t), so that {tilde over (h)} ₂ (t) contains mostlycomponents from the excess channel h ₂ (t). The second constraint inequation (12b) suppresses time variation in the main channel h ₁(t). Thethird constraint in equation (12c) provides an unbiased estimate of h₂(t).

As a specific example, a 3-tap time-domain filter may be used to derivethe 2P channel taps in {tilde over (h)} _(2P×1)(t) based on ĥ(t−1),ĥ(t), and ĥ(t+1) for three symbol periods. The 3-tap time-filter may bedesigned as follows. Using equation (10), the l-th channel tap in symbolperiods t−1, t, and t+1, prior to the time-domain filtering, may beexpressed as:ĥ _(l)(t−1)=h _(1,l)(t−1)+h _(2,l)(t−1)·W _(M) ^(m) ^(t−1) +n _(l)(t−1),ĥ _(l)(t)=h _(1,l)(t)+h _(2,l)(t)·W _(M) ^(m) ^(t−1) +n _(l)(t), for l=1. . . P,ĥ _(l)(t−1)=h _(1,l)(t−1)+h _(2,l)(t−1)·W _(M) ^(m) ^(t−1) +n_(l)(t−1),  Eq (13)

-   -   where ĥ_(l)(t), h_(1,l)(t), h_(2,l)(t), and n_(l)(t) are the        l-th element of ĥ(t), h ₁(t), h ₂(t), and n(t), respectively;        and    -   m_(t−1), m_(t), and m_(t+1) are the interlaces used for pilot        transmission in symbol periods t−1, t, and t+1, respectively.

For the 3-tap time-domain filter for staggering pattern 410 shown inFIG. 4B, with M=8, m_(t−1)=m_(t)−3, and m_(t+1)=m_(t)+3, the constraintsused to select the coefficients for the main channel estimate {tildeover (h)} ₁(t) may be expressed as:

Cancel excess channel: α₁(−1) · e^(−j3π/4) + α₁(0) + α₁(1) · e^(j3π/4) =0, Suppress time variation: α₁(−1) − α₁(1) = 0, and Provide unbiasedestimate: α₁(−1) + α₁(0) + α₁(1) = 1.The first equation above (to cancel the excess channel) is from equation(11a) and has the form: α₁(−1)·W₈ ^(m) ^(t) ⁺³+α₁(0)·W₈ ^(m) ^(t)+α₁(1)·W₈ ^(m) ^(t) ⁻³=0, which may be simplified as: α₁(−1)·W₈³+α₁(0)+α₁(1)·W₈ ⁻³=0, where W₈ ³=e^(−j3π/4) and W₈ ⁻³=e^(+j3π/4).

The solution to the above set of equations for the main channel is givenas:

$\begin{matrix}{\begin{bmatrix}{\alpha_{1}\left( {- 1} \right)} \\{\alpha_{1}(0)} \\{\alpha_{1}(1)}\end{bmatrix} = {\begin{bmatrix}{1 - {1/\sqrt{2}}} \\{\sqrt{2} - 1} \\{1 - {1/\sqrt{2}}}\end{bmatrix}.}} & {{Eq}\mspace{14mu}(14)}\end{matrix}$Equation (14) indicates that the coefficients for the main channelestimate {tilde over (h)} ₁(t) are independent of symbol period t. Thisset of coefficients suppresses time variation in the main channel h ₁(t)but does not suppress time-variation in the excess channel h ₂ (t).Time-variation error is proportional to the energy of the channel taps,which is typically small for the excess channel and significant onlywhen the transmitting and/or receiving entity is moving at high speeds.Thus, not suppressing time variation in the excess channel h₂(t) mayonly marginally degrade performance, if at all.

For the 3-tap time-domain filter for staggering pattern 410 shown inFIG. 4B, the constraints used to select the coefficients for the excesschannel estimate {tilde over (h)} ₂ (t) may be expressed as:

Cancel main channel: α₂(−1) + α₂(0) + α₂(1) = 0, Suppress timevariation: α₂(−1) − α₂(1) = 0, and Provide unbiased estimate: α₂(−1) ·e^(−j3π/4) + α₂(0) + α₂(1) · e^(j3π/4) = e^(j2π·(m) ^(t) ^(−1)/8).The third equation above (to provide an unbiased estimate) is fromequation (12c) and has the form: α₂(−1)·W₈ ^(m) ^(t) ⁺³+α₂(0)·W₈ ^(m)^(t) +α₂(1)·W₈ ^(m) ^(t) ⁻³=1, which may be simplified as: α₁(−1)·W₈³+α₁(0)+α₁(1)·W₈ ⁻³=W₈ ^(m) ^(t) , where m_(t)ε{1 . . . M}.

The solution to the above set of equations for the excess channel isgiven by:

$\begin{matrix}{\begin{bmatrix}{\alpha_{2}\left( {- 1} \right)} \\{\alpha_{2}(0)} \\{\alpha_{2}(1)}\end{bmatrix} = {\begin{bmatrix}{{- 1} + {1/\sqrt{2}}} \\{2 - \sqrt{2}} \\{{- 1} + {1/\sqrt{2}}}\end{bmatrix} \cdot {{\mathbb{e}}^{j\; 2\;{\pi \cdot {{({m_{t} - 1})}/8}}}.}}} & {{Eq}\mspace{14mu}(15)}\end{matrix}$Equation (15) indicates that the coefficients for the excess channel aredependent on the staggering phase m_(t) of interlace m_(t) used forpilot transmission in symbol period t.2. Channel Impulse Response Estimate of Length 3P

To obtain a longer impulse response estimate {tilde over (h)} _(3P×1)(t)with L=3P channel taps, the initial impulse response estimate ĥ(t)obtained in symbol period t for one interlace may be expressed as:ĥ (t)= h ₁(t)+ h ₂(t)·W _(M) ^(m) ^(t) +h ₃(t)·W _(M) ^(2m) ^(t) +n(t),  Eq (16)Equation (16) is derived based on equation (7) and assumes that segments4 through M contain channel taps with zero magnitude. The vectors h₁(t), h ₂(t), and h ₃(t) contain P channel taps for the first, second,and third segments, respectively, of h _(N×1)(t).

A 3-tap time-domain filter may also be used to derive the 3P elements of{tilde over (h)} _(3P×1)(t) based on ĥ(t−1), ĥ(t), and ĥ(t+1) obtainedin three symbol periods. Using equation (16), the l-th channel tap insymbol periods t−1, t, and t+1, prior to the time-domain filtering, maybe expressed in matrix form, as follows:

$\begin{matrix}{\left\lbrack \begin{matrix}{{\hat{h}}_{l}\left( {t - 1} \right)} \\{{\hat{h}}_{l}(t)} \\{{\hat{h}}_{l}\left( {t + 1} \right)}\end{matrix} \right\rbrack = {{\left\lbrack \begin{matrix}1 & W_{M}^{{- \Delta}\; m} & W_{M}^{{- 2}\;\Delta\; m} \\1 & 1 & 1 \\1 & W_{M}^{\Delta\; m} & W_{M}^{2\;\Delta\; m}\end{matrix} \right\rbrack \cdot \left\lbrack \begin{matrix}{h_{1,l}(t)} \\{{h_{2,l}(t)} \cdot W_{M}^{m_{t}}} \\{{h_{3,l}(t)} \cdot W_{M}^{2m_{t}}}\end{matrix} \right\rbrack} + {\quad{\left\lbrack \begin{matrix}{n_{l}\left( {t - 1} \right)} \\{n_{l}(t)} \\{n_{l}\left( {t + 1} \right)}\end{matrix} \right\rbrack,\mspace{79mu}{{{for}\mspace{14mu} l} = {1\mspace{14mu}\ldots\mspace{14mu} P}},}}}} & {{Eq}\mspace{14mu}(17)}\end{matrix}$where

$W_{M}^{\Delta\; m} = {{\mathbb{e}}^{{- j}\; 2\pi\frac{\Delta\; m}{M}}.}$Equation (17) assumes that m_(t−1)=m_(t)−Δm and m_(t+1)=m_(t)+Δm. The3-tap time-domain filter does not have enough degrees of freedom tosuppress time variation in h ₁(t), h ₂ (t), or h ₃ (t). Thus, equation(17) further assumes that h ₁(t), h ₂ (t), and h ₃(t) are constant overthe three symbol periods t−1, t, and t+1.

A least-squares estimate of h ₁(t), h ₂ (t), and h ₃ (t) may be obtainedas follows:

$\begin{matrix}{{\begin{bmatrix}{h_{1,l}(t)} \\{{h_{2,l}(t)} \cdot W_{M}^{m_{t}}} \\{{h_{3,l}(t)} \cdot W_{M}^{2m_{t}}}\end{bmatrix} = {\begin{bmatrix}1 & W_{M}^{{- \Delta}\; m} & W_{M}^{{- 2}\;\Delta\; m} \\1 & 1 & 1 \\1 & W_{M}^{\Delta\; m} & W_{M}^{2\;\Delta\; m}\end{bmatrix}^{- 1} \cdot \begin{bmatrix}{{\hat{h}}_{l}\left( {t - 1} \right)} \\{{\hat{h}}_{l}(t)} \\{{\hat{h}}_{l}\left( {t + 1} \right)}\end{bmatrix}}},{{{for}\mspace{14mu} l} = {1\mspace{14mu}\ldots\mspace{14mu}{P.}}}} & {{Eq}\mspace{14mu}(18)}\end{matrix}$

The 3-tap time-domain filter for h ₁(t), h ₂(t), and h ₃ (t) may beexpressed in matrix form, as follows:

$\begin{matrix}{{\begin{bmatrix}{{\overset{\sim}{h}}_{1,l}(t)} \\{{\overset{\sim}{h}}_{2,l}(t)} \\{{\overset{\sim}{h}}_{3,l}(t)}\end{bmatrix} = {\begin{bmatrix}{\alpha_{1}(1)} & {\alpha_{1}(0)} & {\alpha_{1}\left( {- 1} \right)} \\{\alpha_{2}(1)} & {\alpha_{2}(0)} & {\alpha_{2}\left( {- 1} \right)} \\{\alpha_{3}(1)} & {\alpha_{3}(0)} & {\alpha_{3}\left( {- 1} \right)}\end{bmatrix} \cdot \begin{bmatrix}{{\hat{h}}_{l}\left( {t - 1} \right)} \\{{\hat{h}}_{l}(t)} \\{{\hat{h}}_{l}\left( {t + 1} \right)}\end{bmatrix}}},{{{for}\mspace{14mu} l} = {1\mspace{14mu}\ldots\mspace{14mu}{P.}}}} & {{Eq}\mspace{14mu}(19)}\end{matrix}$

For staggering pattern 410 shown in FIG. 4B, with M=8 and Δm=3, thecoefficients for the 3-tap time-domain filter may be derived based onequation (18) and given as:

${\begin{bmatrix}{\alpha_{1}(1)} \\{\alpha_{1}(0)} \\{\alpha_{1}\left( {- 1} \right)}\end{bmatrix} = \begin{bmatrix}{0.3536 - {j\; 0.1464}} \\0.2929 \\{0.3536 + {j\; 0.1464}}\end{bmatrix}},{\begin{bmatrix}{\alpha_{2}(1)} \\{\alpha_{2}(0)} \\{\alpha_{2}\left( {- 1} \right)}\end{bmatrix} = \begin{bmatrix}{\left( {{- 0.2071} - {j\; 0.2071}} \right) \cdot {\mathbb{e}}^{j\frac{\pi}{4}m_{t}}} \\{0.4142 \cdot {\mathbb{e}}^{j\frac{\pi}{4}m_{t}}} \\{\left( {{- 0.2071} + {j\; 0.2071}} \right) \cdot {\mathbb{e}}^{j\frac{\pi}{4}m_{t}}}\end{bmatrix}},{{{and}\begin{bmatrix}{\alpha_{3}(1)} \\{\alpha_{3}(0)} \\{\alpha_{3}\left( {- 1} \right)}\end{bmatrix}} = {\begin{bmatrix}{\left( {{- 0.1464} + {j\; 0.3536}} \right) \cdot {\mathbb{e}}^{j\frac{\pi}{2}m_{t}}} \\{0.2929 \cdot {\mathbb{e}}^{j\frac{\pi}{2}m_{t}}} \\{\left( {{- 0.1464} - {j\; 0.3536}} \right) \cdot {\mathbb{e}}^{j\frac{\pi}{2}m_{t}}}\end{bmatrix}.}}$

The main channel estimate {tilde over (h)} ₁(t) may be obtained byapplying the coefficients α₁(1), α₁(0), and α₁(−1) to ĥ(t−1), ĥ(t), andĥ(t+1), respectively. The excess channel estimate {tilde over (h)} ₂(t)may be obtained by applying the coefficients α₂ (1), α₂ (0), and α₂(−1)to ĥ(t−1), ĥ(t), and ĥ(t+1), respectively. The excess channel estimate{tilde over (h)} ₃(t) may be obtained by applying the coefficientsα₃(1), α₃(0), and α₃(−1) to ĥ(t−1), ĥ(t), and ĥ(t+1), respectively.

The 3-tap time-domain filter does not have sufficient degrees of freedomto apply many of the constraints shown in equation sets (11) and (12).The coefficients for this time-domain filter do not suppress timevariation in the main channel h ₁(t) or the excess channel hh₂(t) and h₃(t). The various constraints described above may be applied by using atime-domain filter with more than three taps.

In general, a different set of coefficients {α_(s)(i)} may be derivedfor the time-domain filter for the impulse response estimate {tilde over(h)}h_(s)(t) for each segment s. The coefficients for each segment s maybe selected based on various constraints such as: canceling the othersegments, suppressing estimation error due to time variation in thechannel, providing an unbiased estimate of h _(s)(t), minimizing thenoise variance in {tilde over (h)} _(s)(t), and so on. The number oftaps for the time-domain filter determines the number of constraintsthat may be applied to the coefficients. Several exemplary 3-taptime-domain filter designs have been described above. Other time-domainfilters may also be designed based on the description above and arewithin the scope of the invention.

In general, a longer impulse response estimate with L channel taps maybe obtained based on pilot symbols received on L different subbands inone or more symbol periods. The pilot may be transmitted on oneinterlace in each symbol period to limit the amount of overhead for thepilot. The pilot may be transmitted on different interlaces withstaggered subbands in different symbol periods. This allows thereceiving entity to obtain a longer impulse response estimate with morethan P channel taps. A full-length impulse response estimate with all Nchannel taps may be obtained if the pilot is transmitted on all Minterlaces using a complete staggering pattern.

The receiving entity may derive a longer impulse response estimate{tilde over (h)} _(L×1)(t) of length L by filtering initial impulseresponse estimates ĥ of length P for a sufficient number of (S or more)different interlaces. If the pilot is transmitted on a differentinterlace in each symbol period, then the time-domain filtering may beperformed over a sufficient number of (S or more) symbol periods toobtain {tilde over (h)} _(L×1)(t). A progressively longer impulseresponse estimate may be obtained by filtering over more symbol periods.Time-domain filtering over fewer symbol periods provides better trackingof changes in the channel, is thus more robust to Doppler effects, andcan provide an impulse response estimate with a shorter length.Time-domain filtering over more symbol periods increases error in {tildeover (h)} _(L×1)(t) due to changes in the channel over time, is lessrobust to Doppler effects, but can provide an impulse response estimatewith a longer length.

A longer impulse response estimate contains excess channel taps. Sinceeach channel tap contains the complex channel gain at that tap positionas well as noise, a progressively longer impulse response estimatecontains more information regarding the channel but also contains morenoise. The noise from the excess channel taps may be viewed as noiseenhancement resulting from extending the length of the channel estimatebeyond P. If the excess channel energy is relatively small, or if theexcess channel taps are not needed, then better performance may beachieved with a shorter impulse response estimate (e.g., {tilde over(h)} _(2P×1)(t)). If the excess channel energy is relatively large, orif the excess channel taps are pertinent, then a longer impulse responseestimate (e.g., {tilde over (h)} _(3P×1)(t)) may provide betterperformance even with the noise enhancement. Channel estimates withdifferent lengths may be derived and used for different purposes at thereceiving entity.

3. Data Detection

For data detection, a longer impulse response estimate {tilde over (h)}_(2P×1)(t) with 2P channel taps may provide a good trade-off between alonger channel estimate and additional noise from the excess channel.The longer channel estimate mitigates the deleterious aliasing effectshown in equation (7) due to undersampling the frequency domain,provides a more accurate estimate of the main channel h ₁(t), and allowsfor estimation of the excess channel h ₂(t). The longer impulse responseestimate {tilde over (h)} _(2P×1)(t) may be derived as described above.

FIG. 5 shows a flow diagram of a process 500 for deriving a channelestimate used for data detection and decoding. Received pilot symbolsare obtained for the subbands in interlace m_(t) used for pilottransmission in the current symbol period t (block 512). An initialfrequency response estimate Ĥ(t) is derived based on the received pilotsymbols, as shown in equation (5) (block 514). An initial impulseresponse estimate ĥ(t) is derived based on the initial frequencyresponse estimate Ĥ(t), as shown in equation (6) (block 516). Initialimpulse response estimates for at least S₁ symbol periods are filteredwith a time-domain filter having at least S₁ taps to obtain a longerimpulse response estimate {tilde over (h)} _(L) ₁ _(×1)(t) with L₁channel taps, where L₁=S₁·P (block 518).

Post-processing may be performed on the L₁ channel taps in {tilde over(h)} _(L) ₁ _(×1)(t) to further improve channel estimation performance(block 520). The post-processing may include truncation, e.g., settingchannel taps P+1 through L₁ for the excess channel estimate to zeros.The post-processing may alternatively or additionally includethresholding, e.g., setting channel taps in the main and/or excesschannel estimates having energy below a given threshold to zeros. Theunprocessed or post-processed longer impulse response estimate {tildeover (h)} _(L) ₁ _(×1)(t) may then be extended to length N byzero-padding to obtain a vector {tilde over (h)} _(N×1)(t) of length N(also block 520). An N-point FFT may then be performed on {tilde over(h)} _(N×1)(t) to obtain a frequency response estimate {tilde over (H)}_(N×1)(t) for all N subbands (block 522), as follows:{tilde over (H)} _(N×1)(t)= W _(N×N) ·{tilde over (h)} _(N×1)(t).  Eq(20)Process 500 may be performed for each symbol period with pilottransmission.

{tilde over (H)} _(N×1)(t) contains N channel gains for the N totalsubbands and may be expressed as: {tilde over (H)} _(N×1)(t)=[{tildeover (H)} ₁ ^(T)(t){tilde over (H)} ₂ ^(T)(t) . . . {tilde over (H)}_(M) ^(T)(t)]^(T), where {tilde over (H)} _(m)(t) contains P channelgain estimates for P subbands in interlace m. The M frequency responseestimates {tilde over (H)} _(m)(t) for the M interlaces may havedifferent noise variances depending on the particular staggering patternused for pilot transmission. In general, a staggering pattern that ismore spread out (e.g., staggering pattern 410) may result in less noisevariation across {tilde over (H)} _(m)(t) for the M interlaces than astaggering pattern that is more closely spaced (e.g., staggering pattern400).

4. Time Tracking

The receiving entity performs time tracking to estimate and track symboltiming across different OFDM symbols. The symbol timing is used tocapture a window of N input samples (often called an FFT window) fromamong the N+C input samples for each received OFDM symbol. Accuratesymbol timing is pertinent since performance of both channel estimationand data detection is affected by the placement of the FFT window. Thetiming of the received OFDM symbol for each symbol period may beestimated by deriving a longer impulse response estimate for that symbolperiod and detecting for the timing based on an appropriate criterion,e.g. maximizing the energy that falls within the cyclic prefix.

If pilot symbols are available on L different subbands and a timingreference is not available, then a longer impulse response estimate withL channel taps may be derived but only L/2 channel taps may be resolvedwithout any ambiguity. This is because a negative timing error resultsin earlier channel taps aliasing and appearing at the end of the impulseresponse estimate. Thus, it is not possible to determine whether thechannel taps at the end of the impulse response estimate are laterchannel taps (if the symbol timing is correct) or earlier channel tapsthat have aliased (if there is a negative timing error). A longerchannel impulse response estimate with up to N channel taps may beobtained by filtering the initial impulse response estimates for Mdifferent interlaces. The resolvable length of the communication channelis increased by the use of the longer impulse response estimate.

FIGS. 6A and 6B illustrate ambiguity in a channel impulse responseestimate due to timing uncertainty. FIG. 6A shows a channel impulseresponse estimate 610 of length L for an actual channel with an impulseresponse of length greater than L/2. In FIG. 6A, the symbol timing iscorrect and channel impulse response estimate 610 properly includesresponses 612 and 614 of the actual channel at the proper locations.

FIG. 6B shows an impulse response 620 of length greater than L/2 foranother actual channel. If there is no timing error, then a channelimpulse response estimate for this channel would include responses 622and 624 at the locations as shown in FIG. 6B. However, if there is atiming error of x, then response 622 would alias and appear as response632. Thus, a channel impulse response estimate for this channel, withtiming error of x, would be similar to channel impulse response estimate610 in FIG. 6A

FIGS. 6A and 6B illustrate that channel impulse response estimate 610may be obtained for (1) a channel having the impulse response shown inFIG. 6A, with no timing error, or (2) a channel having the impulseresponse shown in FIG. 6B, with a timing error of x, and these two casescannot be distinguished. However, this ambiguity problem would not occurif the channel response lengths are always assumed to be less than L/2.Since the actual channel in FIG. 6B would then have to be longer thanL/2 to be mistaken with the channel in FIG. 6A, it can be concluded thatthe channel response estimate in FIG. 6A does correspond to the truechannel. Thus, an initial estimate of length L can resolve a channel oflength L/2 with timing uncertainty. A longer channel impulse responseestimate is thus desirable for time tracking.

The longer channel impulse response estimate has additional noise due tothe excess channel taps, and greater error due to channeltime-variations. However, time tracking is likely to be less sensitiveto the additional noise since the goal of time tracking is to determineless detailed information such as the general location of the channelenergy rather than the complex channel gains of each tap. Thus, thetradeoff between channel quality and length is consistent with therequirements for data detection and time tracking. Specifically, fortime tracking, a longer impulse response estimate {tilde over (h)}_(3P×1)(t) with 3P channel taps may provide a good trade-off betweenresolvable channel length and noise enhancement. For example, if P=512,then {tilde over (h)} _(3P×1)(t) contains 1536 channel taps, and up to768 channel taps may be resolved without ambiguity. Once the symboltiming is known, the communication channel may be assumed to be 3P/2taps long for data detection purpose. A 3P/2-tap channel may beestimated by obtaining a longer impulse response estimate with 2Pchannel taps and truncating the last 256 channel taps.

FIG. 7 shows a flow diagram of a process 700 for performing timetracking. Blocks 712, 714, 716, and 718 in FIG. 7 are as described abovefor block 512, 514, 516, and 518, respectively, in FIG. 5. However, alonger impulse response estimate {tilde over (h)} _(L) ₂ _(×1)(t) with adifferent length L₂ may be used for time tracking, and a differenttime-domain filter with at least S₂ channel taps may be used to derive{tilde over (h)} _(L) ₂ _(×1)(t), where L₂=S₂·P. The channel estimate{tilde over (h)} _(L) ₂ _(×1)(t) is then processed to determine thetiming of the received OFDM symbol for the current symbol period t(block 720). One method to determine the timing is as follows. A windowof length L₂/2 is placed such that the left edge of the window isinitially at tap index 1. The energy of all channel taps falling withinthe window is computed. The window is then moved to the right, one tapposition at a time until tap index L₂/2 is reached. The channel tapenergy is computed for each tap position. The peak energy among all ofthe L₂/2 window starting positions is then determined. If multiplewindow starting positions have the same peak energy, then the leftmostwindow starting position with the peak energy is identified. Theleftmost window starting position with the peak energy uniquelydetermines the FFT window for the received OFDM symbol. Timing detectionmay also be performed using other techniques. In any case, symbol timingestimate is updated with the timing information obtained for the currentreceived OFDM symbol (block 722).

In general, the same or different impulse response estimates may be usedfor data detection/decoding and time tracking. The use of the sameimpulse response estimate can reduce the amount of computation at thereceiving entity. In this case, the channel length L and the time-domainfilter for this impulse response estimate may be selected to providegood performance for both data detection and time tracking. Differentimpulse response estimates may also be used for data detection/decodingand time tracking in order to achieve better performance for both, andmay be derived with two time-domain filters. The channel length and thetime-domain filter coefficients for each impulse response estimate maybe selected to provide good performance for data detection or timetracking.

FIG. 8 shows an embodiment of OFDM demodulator 160, channel estimator172, and time tracking unit 162 at receiving entity 150. Within OFDMdemodulator 160, a cyclic prefix removal unit 812 captures N inputsamples for each received OFDM symbol based on the symbol timingprovided by time tracking unit 162. An FFT unit 814 performs an N-pointFFT on each window of N input samples and obtains N received symbols forthe N subbands. FFT unit 814 provides received data symbols to detector170 and received pilot symbols to channel estimator 172. Detector 170also receives the frequency response estimate {tilde over (H)} _(N×1)(t)from channel estimator 172, performs data detection on the received datasymbols, and provides detected data symbols.

Within channel estimator 172, a pilot detector 822 removes themodulation on the received pilot symbols and may perform extrapolationand/or interpolation to obtain the initial frequency response estimateĤ(t) composed of P channel gains for the P subbands in the interlaceused for pilot transmission in the current symbol period t. An IFFT unit824 performs a P-point IFFT on Ĥ(t) to obtain the modulated impulseresponse estimate ĥ _(m)(t) with P channel taps. A rotator 826 removesthe phase ramp in the P elements of ĥ _(m)(t) and provides the initialimpulse response estimate ĥ(t). A time-domain filter 830 filters theinitial impulse response estimates ĥ(t) obtained for S₁ or moreinterlaces obtained in S₁ or more symbol periods and provides the longerimpulse response estimate {tilde over (h)} _(L×1)(t) with L₁ channeltaps. A post-processor 832 performs post-processing (e.g., truncation,thresholding, and so on) and zero-padding on {tilde over (h)} _(L) ₁_(×1)(t) and provides a vector {tilde over (h)} _(N×1)(t) with N channeltaps. An FFT unit 834 performs an N-point FFT on {tilde over (h)} _(N×1)(t) to obtain the frequency response estimate {tilde over (H)} _(N×1)(t)for the N total subbands. Channel estimator 172 may also derive afrequency response estimate {tilde over (H)} _(m)(t) for just one ormore selected interlaces.

Within time tracking unit 162, a time-domain filter 840 filters theinitial impulse response estimates ĥ(t) for S₂ or more interlacesobtained in S₂ or more symbol periods and provides the longer impulseresponse estimate {tilde over (h)} _(L) ₂ _(×1)(t) with L₂ channel taps.A timing detector 842 determines the timing for the current receivedOFDM symbol, e.g., based on the energy of the channel taps in {tildeover (h)} _(L) ₂ _(×1)(t). A time tracking loop 844 (which may be a loopfilter) adjusts the symbol timing from the timing used for the currentreceived OFDM symbol.

FIG. 9 shows a block diagram of a time-domain filter 830 x, which may beused for filters 830 and 840 in FIG. 8. Within filter 830 x, the l-thchannel tap in ĥ(t) is provided to N_(f)+N_(b)−1 series-coupled delayelements 912. Each delay element 912 delays its input channel tap by onesymbol period. N_(f)+N_(b)−1 multipliers 914 couple to the input of theN_(f)+N_(b)−1 delay elements, and one multiplier 914 couples to theoutput of the last delay element. The N_(f)+N_(b) multipliers receiveand multiply their channel taps ĥ_(l)(t+N_(f)) through ĥ_(l)(t−N_(b)+1)with coefficients α_(s,l)(−N_(f)) through α_(s,l)(N_(b)−1),respectively. The same coefficients may be used for all P channel tapsin each segment, in which case the coefficients may be denoted asα₈(−N_(f)) through α_(s)(N_(b)−1), without subscript l for tap index. Asummer 916 receives and sums the outputs of all N_(f)+N_(b) multipliersand provides the l-th channel tap in segment s of {tilde over (h)}_(L×1)(t). L may be equal to L₁ for data detection and to L₂ for timetracking. The filtering for only one channel tap in {tilde over (h)}_(L×1)(t) is shown in FIG. 9. The filtering for each of the remainingchannel taps in {tilde over (h)} _(L×1)(t) may be performed in similarmanner.

The pilot transmission, channel estimation, and time tracking techniquesdescribed herein may be implemented by various means. For example, thesetechniques may be implemented in hardware, software, or a combinationthereof. For a hardware implementation, the processing units used forpilot transmission at the transmitting entity may be implemented withinone or more application specific integrated circuits (ASICs), digitalsignal processors (DSPs), digital signal processing devices (DSPDs),programmable logic devices (PLDs), field programmable gate arrays(FPGAs), processors, controllers, micro-controllers, microprocessors,other electronic units designed to perform the functions describedherein, or a combination thereof. The processing units used for channelestimation and time tracking at the receiving entity may also beimplemented within one or more ASICs, DSPs, and so on.

For a software implementation, these techniques may be implemented withmodules (e.g., procedures, functions, and so on) that perform thefunctions described herein. The software codes may be stored in a memoryunit (e.g., memory unit 142 or 192 in FIG. 1) and executed by aprocessor (e.g., controller 140 or 190). The memory unit may beimplemented within the processor or external to the processor.

Headings are included herein for reference and to aid in locatingcertain sections. These headings are not intended to limit the scope ofthe concepts described therein under, and these concepts may haveapplicability in other sections throughout the entire specification.

The previous description of the disclosed embodiments is provided toenable any person skilled in the art to make or use the presentinvention. Various modifications to these embodiments will be readilyapparent to those skilled in the art, and the generic principles definedherein may be applied to other embodiments without departing from thespirit or scope of the invention. Thus, the present invention is notintended to be limited to the embodiments shown herein but is to beaccorded the widest scope consistent with the principles and novelfeatures disclosed herein.

1. A method of performing channel estimation and time tracking in amulti-carrier communication system, comprising: receiving, at areceiver, one or more symbols over a plurality of frequency subbands;obtaining a plurality of initial impulse response estimates for aplurality of different groups of the plurality of frequency subbands,wherein each group includes a different plurality of frequency subbandsselected from among N frequency subbands in the system, where N is aninteger greater than one; processing a first set of initial impulseresponse estimates, from among the plurality of initial impulse responseestimates, to obtain a first impulse response estimate of length L₁ andused for data detection and decoding, where L₁ is an integer one orgreater; and processing a second set of initial impulse responseestimates, from among the plurality of initial impulse responseestimates, to obtain a second impulse response estimate of length L₂ andused for time tracking, wherein L₂ is longer than L₁.
 2. The method ofclaim 1, wherein each of the plurality of groups includes P frequencysubbands, where P is an integer greater than one.
 3. The method of claim2, wherein the P frequency subbands in each group are uniformlydistributed across the N frequency subbands.
 4. The method of claim 2,wherein L₁ is equal to S₁ times P and L₂ is equal to S₂ times P, whereS₁ and S₂ are integers greater than one and S₂ is greater than or equalto S₁.
 5. The method of claim 2, wherein L₁ is equal to two times P andL₂ is equal to three times P.
 6. The method of claim 4, wherein thefirst impulse response estimate comprises L₁ channel taps arranged in S₁segments, each segment including P consecutive channel taps, and whereinthe P channel taps in each segment are obtained based on a set of filtercoefficients for the segment.
 7. The method of claim 6, wherein a firstsegment of P channel taps for the first impulse response estimate is fora main channel, and wherein remaining ones of the S₁ segments for thefirst impulse response estimate are for an excess channel.
 8. The methodof claim 7, wherein the set of filter coefficients for the first segmentof P channel taps is selected to cancel contribution from the excesschannel, suppress time variation in the main channel, provide anunbiased estimate of the main channel, minimize a norm of the filtercoefficients, or a combination thereof.
 9. The method of claim 7,wherein the set of filter coefficients for a second segment of P channeltaps, for one of the remaining ones of the S₁ segments for the firstimpulse response estimate, is selected to cancel contribution from themain channel, suppress time variation in the main channel, suppress timevariation in the excess channel, provide an unbiased estimate of thesecond segment of P channel taps, minimize a norm of the filtercoefficients, or a combination thereof.
 10. The method of claim 4,wherein the first set includes S₁ initial impulse response estimates,where S₁ is an integer greater than one, and wherein the first impulseresponse estimate is derived with a filter having S₁ taps for the S₁initial impulse response estimates.
 11. The method of claim 10, whereinS₁ is equal to two.
 12. The method of claim 4, wherein the secondimpulse response estimate comprises L₂ channel taps arranged in S₂segments, each segment including P consecutive channel taps, and whereinthe P channel taps in each segment are obtained based on a set of filtercoefficients for the segment.
 13. The method of 6, wherein the secondimpulse response estimate comprises L₂ channel taps arranged in S₂segments, each segment including P consecutive channel taps, wherein theP channel taps in each segment are obtained based on a set of filtercoefficients for the segment, and wherein the filter coefficients usedfor the first impulse response estimate are different from the filtercoefficients used for the second impulse response estimate.
 14. Themethod of claim 1, wherein the first impulse response estimate comprisesL₁ channel taps, and wherein the second impulse response estimatecomprises the L₁ channel taps of the first channel impulse responseestimate and (L₂−L₁) additional channel taps.
 15. The method of claim 1,wherein the second set includes S₂ initial impulse response estimates,where S₂ is an integer greater than one, and wherein the second impulseresponse estimate is derived with a filter having S₂ taps for the S₂initial impulse response estimates.
 16. The method of claim 1, furthercomprising: processing the first impulse response estimate to obtain afrequency response estimate for at least one group of frequencysubbands.
 17. The method of claim 16, wherein the processing the firstimpulse response estimate comprises setting channel taps in the firstimpulse response estimate with energy less than a predeterminedthreshold to zero.
 18. The method of claim 16, wherein the processingthe first impulse response estimate comprises setting a predeterminednumber of latest channel taps in the first impulse response estimate tozero.
 19. The method of claim 1, further comprising: determining a timeinstant with a predominant amount of energy in the second impulseresponse estimate; and deriving symbol timing with the time instantdetermined for the second impulse response estimate.
 20. The method ofclaim 1, further comprising: obtaining a plurality of initial frequencyresponse estimates for the plurality of different groups of theplurality of frequency subbands; and processing each of the plurality ofinitial frequency response estimates to obtain a corresponding one ofthe plurality of initial impulse response estimates.
 21. The method ofclaim 20, wherein the processing each of the plurality of initialfrequency response estimates comprises transforming the initialfrequency response estimate to obtain a modulated impulse responseestimate, and rotating channel taps in the modulated impulse responseestimate to obtain the corresponding initial impulse response estimate.22. An apparatus in a multi-carrier communication system, comprising: aprocessing unit operative to obtain a plurality of initial impulseresponse estimates for a plurality of different groups of frequencysubbands, wherein each group includes a different plurality of frequencysubbands selected from among N frequency subbands in the system, where Nis an integer greater than one; a first filter operative to process afirst set of initial impulse response estimates, from among theplurality of initial impulse response estimates, to obtain a firstimpulse response estimate of length L₁ and used for data detection anddecoding, where L₁ is an integer one or greater; and a second filteroperative to process a second set of initial impulse response estimates,from among the plurality of initial impulse response estimates, toobtain a second impulse response estimate of length L₂ and used for timetracking, wherein L₂ is longer than L₁.
 23. The apparatus of claim 22,wherein each of the first and second impulse response estimatescomprises multiple segments, each segment including P consecutivechannel taps, and wherein the P channel taps in each segment are derivedbased on a set of filter coefficients for the segment.
 24. The apparatusof claim 22, wherein the set of filter coefficients for each segment ofeach of the first and second impulse response estimates is selected tocancel contribution from remaining ones of the multiple segments,suppress time variation in the segment, suppress time variation in oneor more remaining segments, provide an unbiased estimate of the segment,minimize a norm of the filter coefficients, or a combination thereof.25. The apparatus of claim 22, further comprising: a transform unitoperative to process the first impulse response estimate to obtain afrequency response estimate for at least one group of frequencysubbands.
 26. The apparatus of claim 22, further comprising: a timingdetector operative to determine a time instant with a predominant amountof energy in the second impulse response estimate; and a time trackingloop operative to derive symbol timing with the time instant determinedfor the second impulse response estimate.
 27. An apparatus in amulti-carrier communication system, comprising: means for obtaining aplurality of initial impulse response estimates for a plurality ofdifferent groups of frequency subbands, wherein each group includes adifferent plurality of frequency subbands selected from among Nfrequency subbands in the system, where N is an integer greater thanone; means for processing a first set of initial impulse responseestimates, from among the plurality of initial impulse responseestimates, to obtain a first impulse response estimate of length L₁ andused data detection and decoding, where L₁ is an integer one or greater;and means for processing a second set of initial impulse responseestimates, from among the plurality of initial impulse responseestimates, to obtain a second impulse response estimate of length L₂ andused for time tracking, wherein L₂ is longer than L₁.
 28. The apparatusof claim 27, wherein each of the first and second impulse responseestimates comprises multiple segments, each segment including Pconsecutive channel taps, and wherein the P channel taps in each segmentare derived based on a set of filter coefficients for the segment.
 29. Acomputer program product, comprising: a computer-readable medium,comprising: code for causing at least one computer to obtain a pluralityof initial impulse response estimates for a plurality of differentgroups of frequency subbands, wherein each group includes a differentplurality of frequency subbands selected from among N frequency subbandsin the system, where N is an integer greater than one; code for causingthe at least one computer to process a first set of initial impulseresponse estimates, from among the plurality of initial impulse responseestimates, to obtain a first impulse response estimate of length L₁ andused data detection and decoding, where L₁ is an integer one or greater;and code for causing the at least one computer to process a second setof initial impulse response estimates, from among the plurality ofinitial impulse response estimates, to obtain a second impulse responseestimate of length L₂ and used for time tracking, wherein L₂ is longerthan L₁.
 30. The computer program product of claim 29, wherein each ofthe first and second impulse response estimates comprises multiplesegments, each segment including P consecutive channel taps, and whereinthe P channel taps in each segment are derived based on a set of filtercoefficients for the segment.